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Recently, several groups have reported observations of collective modes of the charge order present 
in underdoped cuprates. Motivated by these experiments, we study theoretically the oscillations of 
the order parameters, both in the case of pure charge order, and for charge order coexisting with 
superconductivity. Using a hot-spot approximation we find in the coexistence regime two Higgs 
modes arising from hybridization of the amplitude oscillations of the different order parameters. 
One of them has a minimum frequency that is within the single particle energy gap and which is a 
non-monotonic function of temperature. The other - high-frequency - mode is smoothly connected 
to the Higgs mode in the single-order-parameter region, but quickly becomes overdamped in the case 
of coexistence. We explore an unusual low-energy damping channel for the collective modes, which 
relies on the band reconstruction caused by the coexistence of the two orders. For completeness, we 
also consider the damping of the collective modes originating from the nodal quasiparticles. At the 
end we discuss some experimental consequences of our results. 


I. INTRODUCTION 

Despite the decades of intense research efforts, super¬ 
conductivity in cuprates remains a profound mystery. 
However, recently there has been a lot of progress in 
clarifying and refinin g th e phase diagram of these ma¬ 
terials experimentally.1221 In particular, there is growing 
evidence of a charge order existing in the pseudogap state 
of several cuprate families,!^® which also coexists and 
competes with superconductivity at lower temperatures. 
Furthermore, it app ears that this order has a non-trivial 
d-wave phase factor,!^ implying that within a one-band 
model of copper sites it describes ordering entirely on the 
links. For this reason it has been dubbed “bond-density 
wave” (BDW). Recently, several groups employed time- 
domain reflect ivitj®! as a tool to study this order,lii^ti^ 
and, in particular, its collective modes. In some cases 
they were able to extract the amplitude and phase oscil¬ 
lations and to track them as the system became super¬ 
conducting. These results can provide valuable insights 
into the physics of both pseudogap and superconducting 
states, and, thus, it is desirable to have a better the¬ 
oretical understanding of the possible collective modes 
of these systems. One particularly interesting point is 
that the coexistence of charge order and superconductiv¬ 
ity makes possible the direct observation of the supercon¬ 
ducting Higgs mode, as first pointed out in the pioneering 
work of Littlewood and Varma.^ 

In this work we present a theoretical study of the Higgs 
modes, or oscillations of the amplitude, of the order 
parameters in underdoped cuprates. We consider both 
the pure BDW state, as well as the coexistent BDW- 
superconductivity phase. We use the so-called “hot-spot” 
modePH22l of the pseudogap phase, which is based on a 
picture of a metallic state close to a magnetic instabil¬ 


ity, and considers the physics of the special points on the 
Fermi surface connected by the magnetic ordering vector. 
Although relatively simple, this model has seen extensive 
use recently, as it naturally leads to coexistence between 
BDW and superconductivity, and also correctly predicts 
the d-wave phase factor of the charge order.E^Hf^ 

Our results provide a general framework for identify¬ 
ing and understanding order-parameter collective modes 
of the system. In the single-order phase (i.e., only BDW 
or superconductivity) we find, as expected, a single am¬ 
plitude mode, which is coupled with the quasiparticle 
continuum and is always damped. However, the coex¬ 
istence regime is much more interesting - the fluctu¬ 
ations of the different order parameters become inter¬ 
twined.!^ As a consequence, in this region we hnd two 
Higgs modes, which represent coupled oscillations of the 
order parameters. One of the modes is slow, with fre¬ 
quency well below the amplitude of the order parameters, 
but which is, nevertheless, weakly damped. The other 
mode is pushed inside the high-energy continuum, and 
quickly becomes overdamped. We follow the slow mode 
in the entire coexistence phase, and find its frequency 
to be a non-monotonic function of temperature. This 
mode is weakly damped through an unusual low-energy 
decay channel for the antinodal qausiparticles, caused 
by the coexistence of the two orders and the associated 
band reconstruction. Even more unusually, this damp¬ 
ing initially increases with the decrease of temperature. 
To account for the damping from the gapless degrees of 
freedom present at the nodal regions we develop a phe¬ 
nomenological time-dependent Ginzburg-Landau theory. 
We demonstrate that, by allowing for signihcant damp¬ 
ing, the in-gap mode is strongly suppressed, while the 
frequency of the high-energy mode is brought down. Our 
results provide a characterization of the amplitude modes 
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of the coexistent superconductivity-BDW system which 
can be compared with the experimental data and used to 
identify the appropriate Higgs modes of the system. 

II. MICROSCOPIC CALCULATION OF 
COLLECTIVE MODES FREQUENCIES 

We will consider here the collective modes in a 2D 
“hot-spot” model. Such a model can be obtained as a 
low-energy theory from the 2D t-J-V modelE^®^ which 
contains hoppings t(i/ 2 / 3 ) on a square lattice as well as 
nearest-neighbor exchange and Coulomb interactions J 
and V. Specifically, one projects the lattice theory onto 
regions in the vicinity of 8 “hot spots” where the Fermi 
surface intersects the magnetic Brillouin zone boundary. 
In the vicinity of these hot spots the nearest-neighbor 
interactions J and V can be approximated by constants. 
Time reversal symmetry allows the problem to be re¬ 
duced to considering fermions near 4 inequivalent hot 
spots where, in the channels of interest, the interactions 
take the form 

= I E (1) 

k,p,q 

= I E (2) 

k,p,q 

where "^a,b are Nambu spinors in pairs of hot-regions sep¬ 
arated by the antiferromagnetic wave-vector K = (tt, tt) 
and gs and gc are the non-retarded components of the in¬ 
teraction in the superconducting and bond-density wave 
(BDW) channels, respectively. Va and are the vertices 
for pairing in the superconducting and bond-density- 
wave channels. Their explicit forms for the system stud¬ 
ied here are shown in Eq. 

Due to the d-wave symmetry of the order parameters 
one can further restrict attention to 2 of the 8 hot re- 
gions.li^ The interaction terms can be decoupled via a 
Hubbard-Stratonovich transformation. In the usual man¬ 
ner, the saddle point of the zero-frequency terms of the 
decoupling fields leads to a mean-field theory, which in 
this case has mean-field Hamiltonian 

n = + + (3) 

^ 9s 9c 

k 

where now 4^ is a Nambu spinor 
{cki^,Ck 2 ,i,cik 2 i’'^-kit)'^ describing one pair of hot 
spots. The mean-field Hamiltonian describes two species 
(denoted 1 and 2 ) of spinful fermions which pair only 
with each other. Specifically, 

Hmf = Hq + AVa + 

Ido = diag(^i,C2) 0 Tz, (4) 

VA=Pi®h, H0 = /3o(g)fi, 

where Ti and pi are Pauli matrices acting in particle-hole 
space and species space, respectively, and A describes 


d-wave superconductivity while (j) is the BDW order.l^ 
Here, and in what follows, M denotes a matrix in the 
4x4 Nambu-hot-spot space, and M a 2 x 2 matrix. The 
self-consistency equations associated with Eq. are 

A = ^TEtriAGfe, 

^ (5) 

k 


where Gk is the matrix Matsubara Green’s function of the 
Hamiltonian in Eq. and k = {ien, k), with e„ being a 
fermionic Matsubara frequency. Here we have considered 
A and (p to be real and non-negative (they can always be 
brought to this form via a gauge transformation). 

In the case of a hot-spot model of cuprates, the two 
species correspond to fermions within a vicinity of in¬ 
equivalent “hot spots” in the Brillouin zone. Close to 
the hot-spot points the electron dispersion can be mod¬ 
eled as ^i{k) = ^ 2 {—k) = Vfkx + iky, where we include 
the curvature 7 as it plays an important role in break¬ 
ing the deg enera cy between the two orders and allowing 
coexistence! ^^ * ^^ ! 

We follow Ref. [TH] by choosing units where Vf = 1, 
7 = 1/A = TT, with A being the hot spot cutoff, and 
parametrize {go5s} = 3J ± AV with J = 1.2. Note 
that V strengthens the interaction in the charge channel, 
while decreasing the interaction in the superconducting 
channel, and thus can be used to tune the coexistence 
(as depicted in Eig. [^. We consider two qualitatively 
different cases of coexisting charge order and supercon¬ 
ductivity as depicted in Fig. one where charge order 
disappears for some finite temperature below the super¬ 
conducting Tc (V = 0.2), and one where charge order 
survives all the way down to T = 0 (U = 0.21). In both 
of these cases, the BDW order will onset at a temper¬ 
ature Tbdw > Tc- The competition between the two 
orders can be readily confirmed by a decrease in (j) below 
the superconducting T^. 

The ordering vector of the BDW is determined by 
the separati on in th e Brillouin zone of the hot spots 
being paired.R is important to note that we 
are considering a BDW with ordering vector {Q,Q), 
which is k nown to be the leading instability of this sim¬ 
ple model.fi^ElHII] This is different from the experimen¬ 
tally observed bond-oriented ordering directions {Q, 0) 
and (0, Q), which correspond to a different choice of hot 
spots for the BDW pairing to occur betwe en. It is possi¬ 
ble to stabilize the {Q, 0) and (0, Q) orders ,14^1111211 but at 
the price of significantly complicating the model, and we 
will not pursue these modifications here. We expect that 
most of our results and conclusions are applicable to the 
{Q, 0)/(0, Q) orders as well. 
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FIG. 1. (C olor o nline) Schematic phase diagram of the hot¬ 
spot modelp^^ which illustrates the transition from super¬ 
conductivity to charge order, tuned by V (nearest-neighbor 
Coulomb interaction). In this work we consider the transition 
from BDW to BDW-superconducting mixed state along the 
“trajectories” indicated by the dashed lines. Depending on 
the exact value of V the T —> 0 limit of the system could be 
either in a pure superconducting state (indicated by the blue 
dashed line), or in a mixed state (red dot-dashed line). 


A. Hybridized Higgs modes 

The collective modes of coexisting charge-density-wave 
and supercondu cting states have been studied theoreti¬ 
cally previously,fi2l2Mini and we apply the methods de¬ 
veloped in these earlier works. In general, the collective 
modes of the system are described by a 5 x 5 matrix, 
which includes the amplitude and the phase modes of 
each order parameter, as well as the density oscillations 
of the fermions. However, this matrix factorizes into two 
decoupled sectors,!^ with a 2 x 2 block describing the in¬ 
teracting amplitude modes, and the other -3x3- block 
describing the order parameters phases coupled to each 
other, as well as to the fermionic density.!^ This being the 
case, we devote our attention to the amplitude mode sec¬ 
tor. In particular, we consider amplitude fluctuations of 
these order parameters with finite frequency cu, but zero 
wave-vector. Doing so allows us to calculate the mass of 
the collective modes: the minimum energy required to 
excite the collective modes of the ordered state. 

Returning to the Hubbard-Stratonovich decoupling of 
the hot-spot model’s interactions, inclusion of the finite- 
frequency components of the decoupling fields leads to 
the action 

k,q,uJm 

+ X! + IT X! (6) 

ys ^ yc ^ 

q^UJm 

where ^mf is the action corresponding to Eq. and we 
are working in imaginary time. We have kept here the 


fluctuations A(t), (/)(r) which are along the direction of 
A, (p in the complex plane,!^ corresponding to the ampli¬ 
tude modes .1^ 

Particularly we will be interested in the 2x2 matrix 
collective mode propagator 




where 02 ,uj^,q = 4>ujr,,,q, and the 

related object q) = Dijiiujm —t w -I- z0+, g), which 

can be obtained via analytic continuation. The off- 
diagonal elements of this matrix are in general non-zero 
and this is what leads to the hybridization of collective 
modes. The poles of the retarded propagator will 
describe the on-shell collective mode energies. 

After integrating out the fermionic degrees of freedom, 
D{uim,^ can be expressed (at the quadratic level) as 




where we have defined 


Z)° = 


I 

4 


gs 0 
0 Qc 


( 8 ) 

( 9 ) 


Here Qij, also a 2 x 2 matrix, is the self-energy of 
the collective modes due to the fermionic quasiparticles 
(this treatment is equivalent to obtaining the generalized 
susceptibilities of the order parameters within the RPA 
approximation). Since, is already known, Q is the 
object of interest. 

Specifically, Q is given by 


Qij (^^ 


i,g) = -T^tr 

k,en 


G{k, en)V^G{k - g, e„ - uJm)Vj 


m 

where S {A,(^}. After performing the fermionic 
Matsubara sums in Eq. we analytically continue the 
bosonic frequency to the real axis, in order to obtain the 
finite-temperature, retarded self-energy Q^{uj,q). 

The long-wavelength frequencies of the amplitude 
modes are given by the solutions of 


det[(I)°)-i - Q«(cco - iPo, 0)] = 0, (II) 


where 


uiQ = Re[a;(g —>■ 0)], (12) 

Po =-Im[a;(g0)], (13) 

are, respectively, the mass and the decay rate of the Higgs 
mode in the long wavelength limit. The in-gap collective 
modes, are those for which ojq < 2min(((), A). 

One can explicitly show that the diagonal components 
of reproduce the usual 2(j)/2A amplitude modeJ^^in 
the limit where one of the order parameters vanishes. 
However, in our case, we focus our attention on the 
eigenmodes of the response function, which describe hy¬ 
bridized modes of the systenP^ and which cannot be ob¬ 
tained from purely considering the superconducting and 
BDW susceptibilities. 
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Because we are interested in weakly damped oscilla¬ 
tions such that it makes sense to describe them as col¬ 
lective modes, we are able to employ a technique to de¬ 
termine the complex frequency of the oscillations from 
considerations of the response function on the real fre¬ 
quency line. In particular, we obtain the real part of the 
frequency as the solution to the equation Re[A(a;o)] = 0 
where A is a solution to the eigenvalue problem 



Q^{uj) - X{uj)i 



(14) 


The imaginary part of the frequency can then be calcu¬ 
lated by expanding the eigenva lue a s a function of com¬ 
plex w about the real frequency.^^Sl^ We defer analysis of 
the imaginary part (shown in Fig. until Sec. IIB and 
focus now on the real part. 

In order to track the temperature dependence of the 
collective modes, we explicitly solve the mean held equa¬ 
tions for a range of temperatures and then calculate the 
collective mode frequencies at each temperature. Below 
TbdWj in the pure BDW phase, we Hnd an amplitude 
mode starting at frequency 2(j), as expected.^ With the 
onset of superconductivity, another mode appears inside 
the gap. Physically, it represents coupled oscillations of 
the two order parameters, wherein pairs are excited in 
both the BDW and superconducting channels. The mix¬ 
ing of the two orders arises due to the off-diagonal ele¬ 
ments of Q^^ proportional to t/fA. Intuitively, one might 
anticipate the presence of such an in-gap mode by argu¬ 
ing that one could convert one type of pairing into the 
other at a smaller energy than it would take to completely 
break a pair. 

The temperature dependence of the mode’s frequency 
is non-trivial - initially it grows, but then reaches a max¬ 
imum and goes down with the decrease of either ^ or A. 
Depending on the shape of the coexistence region, this 
mode either survives all the way down to T = 0 or it van¬ 
ishes at the second transition to a single-order-parameter 
phase. This behavior can be seen in Fig. Note that 
near the phase transitions, this mode approaches the 
20 / 2 A amplitude mode of the order that vanishes at that 
temperature, which is the reason for the softening of the 
mode in the vicinity of these points. 

At the onset of the coexistent phase, the other (20) 
mode is pushed to higher energies, enters the quasiparti¬ 
cle continuum, and quickly becomes overdamped. Thus, 
it is outside the region of validity of our method of finding 
w, and so we do not track it. 


B. Damping from antinodal quasiparticles 

As explained in Sec. |II A[ the damping rate Fq, can be 
obtained by expanding the eigenvalues of Eq. about 
the real part of the zero momentum dispersion ojq. The 
temperature dependence of this damping rate is shown 
in Fig. 1^ Although the in-gap mode stays below the 
(2A, 20) threshold, its frequency has a finite damping 
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FIG. 2. (Color online) Mass of the in-gap hybrid Higgs mode 
LOO = Re[a ;(5 —>■ 0)] as obtained from Eq. 11 The frequency is 


plotted as a function of temperature for two different cases of 
0(r —>■ 0) (F = 0.2,0.21) as depicted by the dashed lines in 
Fig.B nsing the units of Ref. 1181 A soft mixed mode emerges 
in both cases below the superconducting T^. For reference, 
twice the single-particle energy gap, which is determined by 
2min(A, 0), is plotted in the black dashed line. In proximity 
of a phase transition, the in-gap mode approaches the 2A/20 
Higgs mode of the vanishing order. 


rate, which, furthermore, initially increases as tempera¬ 
ture goes down. This unusual behavior of the damping 
arises from the BDW bubble when just charge or¬ 
der is present, the only scattering which could lead to 
damping requires at least energy 20 (as can be seen in 
Fig. 4a). All other types of scattering have zero ma¬ 
trix element, and thus there is no damping at 9 = 0 for 
Wq < 20. However, as soon as A becomes non-zero the 
bands are reconstructed due to hybridization of the BDW 
bands with their corresponding hole bands, and simul¬ 
taneously scattering matrix elements between all bands 
become non-zero, allowing transitions between any two 
bands to contribute (c.f. Fig. 4b). As a result, there now 
exist transitions for arbitrarily small frequency (between 
the two particle/hole bands), giving rise to the damping 
of collective modes within the gap. 

The specific temperature dependence of the damping 
results from a combination of two effects. Because we 
are considering energies cog < 2min(0, A), we see that 
transitions from a particle to a hole band (or vice versa) 
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FIG. 3. (Color online) Damping rate Fq of the in-gap collec¬ 
tive mode in the long wavelength limit for two different values 
of V corresponding to the two different trajectories depicted 
in Fig. Damping is an order of magnitude smaller in the 
case where <f)(T —>■ 0) 7 ^ 0, and is exponentially suppressed at 
low temperature due a lack of thermally excited quasiparti¬ 
cles. In both cases, the decay rate is strongly suppressed in 
the vicinity of the superconducting T^. The dashed vertical 
line indicates the boundary between the coexistent and pure 
superconductivity phases for the case of the blue curve (c.f. 
the upper plot of Fig. [^. 


cannot contribute as they will always have energy equal 
or greater than 2A. Thus, damping must be solely due 
to scattering between the hole or particle bands. As 4> 
decreases, the two particle (and correspondingly the two 
hole) bands become more similar (in the limit ^ 0 they 

are degenerate), increasing the phase space for low energy 
transitions and therefore leading to greater damping of 
the BDW amplitude mode. This in turn leads to an 
increased damping of the mixed mode, which is visible in 
Fig.i However, in opposition to this effect, as (^ —f 0, 
the matrix element for scattering between these bands 
will begin to vanish, as it is proportional to </>. At some 
point this second effect will overcome the increase due to 
the larger phase space, leading to a disappearance of the 
damping as we approach the critical point at which the 
charge order disappears. 

The competition between scattering elements and band 
structure generically leads to a non-monotonic temper¬ 
ature dependence of the damping, which in turn means 
that there exists a region of maximal damping away from 
which the decay rate remains weak (within the gap). In 
the case with (j)(T —>■ 0) 7 ^ 0, the BDW order remains 
sufficiently large that the system never approaches this 
region of larger damping and thus the decay rate is no¬ 
ticeably smaller than for (j){T —>■ 0) = 0. In all cases 
where the mixed phase exists down to T = 0, this damp¬ 
ing term will be exponentially suppressed at low temper¬ 
atures as there are no thermally excited quasiparticles 
available to scatter. 



K 

(b) 

FIG. 4. (Golor online) Slices of the quasiparticle band struc¬ 
ture for ky = 0. Dashed lines indicate the particle hole 
conjugate of the band of the same color, (a) The normal 
state (gray) and BDW state (green/blue) dispersions. Only 
transitions between the two solid/dashed bands contribute 
to \raQ^^. The onset of superconductivity hybridizes the 
green/blue bands with their particle-hole conjugates, (b) Bo- 
goliubov band dispersions in the coexistent state. Transitions 
between all bands may contribute to ImQ^,^ leading to damp¬ 
ing of the Higgs modes (even for those with mass less than 
2A.). Processes indicated by horizontal arrows are of par¬ 
ticular importance as they occupy a finite phase space for 
arbitrary frequencies within the gap. 


III. DAMPING FROM NODAL 
QUASIPARTICLES 

The hot-spot model we have used so far is only defined 
in the antinodal regions, and thus completely ignores the 
gapless degrees of freedom existing close to the nodes. 
These can have a particularly strong effect on the damp¬ 
ing of the collective modes by providing a low-energy 
decay channel. However, the contribution of these quasi¬ 
particles is different for the different orders. We expect 
the charge order to couple only weakly to the nodal quasi¬ 
particles, due to the mismatch between its wavevector 
(Q, Q) and the wavevector separating the node^ [note 
that the same argument applies to BDW with ((3,0) or 
( 0 ,( 3 ) wavevector]. There is no such restriction for the 
superconductivity, however, and its amplitude fluctua- 
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tions are unavoidably damped by the nodal excitations. 
To include these effects and to study their consequences 
for the collective modes, we supplement the calculation 
from the previous section with a phenomenological time- 
dependent Ginzburg-Landau theory. In addition to the 
more familiar quadratic and quartic in the order param¬ 
eters terms, this theory contains also first and second 
derivatives in (real) time. The time-dependent Ginsburg- 
Landau equations can be written in the following formPZl 


92 A 

~W 


dA dJ=GL 

dA* ' 
dJ^GL 


dt 

d(j) 


(15) 


d(j)* 


where the Ginsburg-Landau action is given by: 


Fgl = a, 


■-l-a AIA p-1-,001 (/) I ^-l-AIAI ^-f u I (/) I ^ IA p. 


The quadratic coefficients a have the usual linear- 
in-temperature dependence, whereas (3 and u (which 
parametrises the competition between the two orders) 
are temperature-independent Note that expressions 
for the coefficients in jFgl can be straightforwardly de¬ 
rived from the microscopic theory presented in the pre¬ 
vious sectiorP^ (spatial derivative terms are not included 
since we are considering only uniform states). Although 
the Ginzburg-Landau theory is strictly applicable only 
close to the critical region, it can be used beyond its 
region of validity as an effective model for the collec¬ 
tive modes of the system.ESl For this reason we keep the 
second-order time derivative terms, which are usually 
omitted close to the critical temperature.^^ 

The coefficients 7 ,^ and 7 a are responsible for the 
damping of the collective modes. It is important to note 
that despite the symmetric way these terms enter Eq. 15 
they encode very different physics. The 7 ^ term is native 
to the hot-spot regions. At low energies it is proportional 
to A, since it is only allowed by the band reconstruc¬ 
tion (see the discussion in the previous section), whereas 
above 2min(A, c^) we can treat it as a constant, origi¬ 
nating from the coupling of the fluctuations to the high- 
energy quasiparticle continuum. In contrast, the main 
contribution to the 7 a term originates from the nodal 
regions (and thus is completely absent in the hot-spot- 
only approach of the previous section). Close to we 
can obtain its temperature dependence from the follow¬ 
ing qualitative considerations. This term is proportional 
to the number of available states at the oscillation fre¬ 
quency, given by ~ p(a;) tanh (a;/4T).l^ Linearizing the 
density of states close to the nodes p{(jj) ~ w, and ap¬ 
proximating the frequency as w « 2A we finally get for 
the damping terms of the slow mode 

7A « 7aA2 = - T), 7 ^ « 7 OA = 4VTc-T 

(we have expanded in powers of A). Note that we have 
thus determined the temperature dependence of 7 ^ and 


7 a, but their relative strength at some fixed tempera¬ 
ture depends on the parameters of the microscopic mod¬ 
els (like V), which cannot be estimated within our phe¬ 
nomenological theory. However, given the general tem¬ 
perature dependence of 7 ^ and 7 a , we expect the antin- 
odal particles to dominate damping sufficiently close to 
Tf. (A vs. A^), whereas at low temperatures the nodal 
excitations take over - 7 a stays finite for T —>■ 0 , while 
70 goes to zero exponentially. 

To obtain the frequencies and damping of the mixed 
modes we expand ^(t) and A(t) around the mean field 
values of the order parameters ^0 and Aq: 4’{t) = 
(j)o + and A{t) = Aq -I- SA{t). Assuming that S(j){t) 
and SA{t) are relatively small we can simplify Eq. 15 by 
keeping only the terms linear in d4> and SA. Since we 
are interested in the collective modes we write their time 
dependence as Inserting this ansatz in the lin¬ 

earized equations, we can exclude S(f> and 6A altogether, 
and finally arrive at the following equation for ui: 


-1-170 w-|-2(Q;0-|-uAg)— 


{2u(l)oAoy 


2(q;a + u4>l) + i7AW -I- w' 


= 0 . 


(16) 


We solve it numerically (with (j)o{T) and Aq{T) deter¬ 
mined by the time-independent mean-field equations), 
and obtain both complex and purely imaginary solutions 
for oj. The former solutions are oscillatory (with Re[a;] 
giving the frequency of the uniform oscillations around 
the mean field values), while the latter represent expo¬ 
nential decay. We show the real and the imaginary parts 
of the two oj solutions as a function of temperature in Fig. 

There we plot wq and Tq for two different strengths 
of 7 ^, as a comparison between small and large contri¬ 
bution from the nodal quasiparticles, respectively. For 
small 7 ^ we can see that both the real and imaginary 
parts of the frequency of the hybridized modes show be¬ 
havior similar to that obtained in the previous section. 
However, when we increase 7 ^ we see not only enhance¬ 
ment of the damping of both modes, but also decrease 
of their real frequencies (the top panel of Fig. [^. Al¬ 
though the effect is more dramatic for the in-gap mode, 
which now exists only in a narrow region below Tc, it is 
significant for the fast one as well. This is a consequence 


of one important feature of Eq. 16 - the coupling of the 
two channels mixes their real and imaginary parts. Thus, 
increase of the damping leads to the gradual suppression 
of the real part of both mixed modes. Note also that the 
disappearance of wq of the in-gap mode corresponds to a 
peak in its Tq. 


IV. DISCUSSION AND CONCLUSION 

Note that our calculation is to some extent complemen¬ 
tary to those in Refs. |23] and |43l These works studied 
the dynamics of the system after an external perturba¬ 
tion, and were done in the time domain, thus allowing 
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T^^bdw 


FIG. 5. (Color online) The evolution of cjq (top panel) and 
To (bottom panel) of the two Higgs modes with temperature. 
For each mode the cases of weak and strong damping from 
the nodal regions are shown [ 7 A = O.la (dashed line) and 
7a = 0.7a (solid line), respectively]. 70 , the damping from 
antinodal region, is the same on both plots. 


direct comparison with the experimental data. The tem¬ 
perature dependence of the frequencies extracted in Ref. 
03] appears consistent with our calculation, as it shows a 
low-frequency mode appearing below the superconduct¬ 
ing transition. 

The experiments have not observed a soft mode close 
to either charge or superconducting transition tempera¬ 
tures. Instead, the frequency of the identified amplitude 
mode stays almost constant, with only a small decrease 
in frequency at the superconducting Tc observed in Ref. 
H and no clear change seen in Ref. [T31 This appears 
consistent with the behavior of the high-energy mode in 
the case of strong damping from the nodal regions (see 
section |Tn|. This damping can effectively “squeeze” the 
low-frequency mode inside a very narrow region close to 
Tc (where it would be difficult to observe), and could also 
lead to the decrease of the frequency of the fast mode, ob¬ 
served in Ref. [12] (note that a different phenomenological 
explanation for this decrease, based on time-independent 
Ginzburg-Landau theory, was given in Ref. mi). In con¬ 
trast, the absence of softening close to Tbdw appears 
incompatible with our calculation, and requires alterna¬ 
tive explanations (such as optical phonons).^ 

In conclusion, we have studied the collective modes 
for the bond density wave and superconducting order 


parameters expected to exist in the pseudogap state of 
cuprates. In the pure BDW phase we observed the con¬ 
ventional amplitude mode with frequency starting at 2(j). 
In the coexistent phase two collective modes represent¬ 
ing the coupled oscillations of the amplitudes of the or¬ 
der parameters are present. One of them is soft at the 
superconducting critical temperature, and despite hav¬ 
ing frequency wq < 2 min(((). A) is (weakly) damped, due 
to band-structure reconstruction caused by superconduc¬ 
tivity. The other mixed mode is continuously connected 
to the pure BDW mode, with frequency pushed up in 
the coexisting regime. To study the effects of damping 
originating from the nodal regions, we developed a phe¬ 
nomenological time-dependent Ginzburg-Landau theory. 
We demonstrated that strong damping can have signifi¬ 
cant effect on the real frequency of the modes. 
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Appendix: Effect of phonons 


Beyond just the non-retarded interaction considered 
above, one can also consider the effect of phonons on 
the collective modes. Here we will take this into account 
by considering the contribution of the frequency depen¬ 
dent phonon-mediated interaction between electrons to 
the collective mode propagators. In particular, we will 
project this interaction onto a hot spot model by taking 
the phonon momentum to be the fixed wavevector Q sep¬ 
arating the hot spots which are being paired - this is the 
same approximation that one uses on the non-retarded 
interaction in deriving the hot spot model. 

A simple Aig symmetry phonon has no effect on the 
collective modes due to the pure d-wave symmetry of the 
order parameters. However, in reality we expect some 
direct order parameter-phonon coupling, either because 
there is a phonon mode with the correct symmetry (Big), 
or because in real systems the order parameter would 
not necessarily have a pure d-wave symmetry, but could 
have an s-wave component admixed. Regardless of the 
exact nature of the coupling, it gives rise to a term in the 
mean field theory which includes the phonon-mediated 
interaction as 


Hph = f Y. 

X (^(l){uJrn)c\„{k, e„ - UJm)c2a{k, £„) -f h.C.^ 


(A.l) 


where 


U {(-Orn) — 


^ep 











T 

FIG. 6. (Color online) The in-gap collective mode mass ujq 
(relative to the minimum pair creation energy) as a function 
of the electron phonon coupling g, for Qq — 1.. The g = 0 
line reproduces the behavior shown in Fig. 

is an Einstein phonon type propagator and / is a constant 
of order one arising from the form factor of the electron- 


phonon vertex. 

If we consider the effect of this term on the charge 
collective mode, we find that it can be captured by the 
replacement gc —>■ gdiuJm)- We absorb the w = 0 compo¬ 
nent into the definition of gc (as that is what determines 
the static mean-field solution) and include the remain¬ 
ing frequency dependent part in our calculation of the 
collective modes. Upon analytic continuation to real fre¬ 
quency, this amounts to the substitution 

gc = 9c- ^ - /gep ^2 ^ ^2 (^-2) 

in the collective mode equations (the additive constant 
is chosen so that we recover gc(w = 0) = gc)- The previ¬ 
ous analysis can now be repeated for a range of electron 
phonon couplings. As can be seen in Fig. the cou¬ 
pling to phonons tends to push the collective mode gap 
downward, while leaving the softening at the phase tran¬ 
sitions unmodified. Overall, the qualitative behavior of 
the mode is not markedly different. 
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